The polymer theta-point as a knot delocalisation transition 
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We study numerically the tightness of prime flat knots in a model of self-attracting polymers with 
excluded volume. We find that these knots are localised in the high temperature swollen regime, but 
become delocalised in the low temperature globular phase. Precisely at the collapse transition, the 
knots are weakly localised. Some of our results can be interpreted in terms of the theory of polymer 
networks, which allows to conjecture exact exponents for the knot length probability distributions. 
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The presence of knots in single ring macromolecules, or 
of chain entanglement in polymer melts, has fundamen- 
tal consequences at the physical, chemical and biological 
level. For example, the dynamics of a ring polymer in a 
gel depends crucially on its topology . The replication 
and the transcription of circular DNA are controlled by 
enzymes affecting the topology, the topoisomerases ||. 
Knots have even been identified in some proteins in their 
native state || . ft is easy to imagine the importance that 
the degree of localisation of knots can have. In the folding 
of a protein it must make quite a difference whether the 
knot is tight and localised within a restricted region of 
the backbone, or not. A loose knot poses less restrictions 
to the exploration of configuration space in the search 
for the native state. Similarly, one could expect that 
the degree of knot localisation can strongly influence the 
function of topoisomerases. 

These examples all concern heteropolymers in nonequi- 
librium situations. Here we investigate the interplay be- 
tween topology and temperature in a simpler context, by 
studying the equilibrium properties of (prime) knots in 
homopolymers when these are cooled below their theta 
temperature. Under such conditions we find that knots 
definitely loose their usual property of being localized 
within restricted portions of the chains. 

It is very difficult to include topological constraints 
within a statistical mechanical description of a polymer, 
since they imply a global control of its conformations j| . 
So far, most work has concentrated on the probability 
of occurence of knots || and much less has been done 
on the physically relevant problem of precisely quanti- 
fying the knot size. Besides attempts to measure knot 
size directly Q, most studies have given only indirect, 
and often incomplete information on knot localisation. 
For a ring polymer in the infinite temperature, good sol- 
vent regime, it was found numerically that the presence 
of a prime knot leads to a simple multiplication of the 
partition sum with a factor proportional to the length 
L of the macromolecule Q]. Moreover, for the relation 
between the radius of gyration and L (Eq.([l])), evidence 
was given that neither the amplitude A |8|] nor the crit- 



ical exponent v depends on the topology [&[|. These 
results suggest that knots are somehow localised within 
the chain. On the other hand, in a recent study of the 
force-extension relation for knotted polymers in good sol- 
vent, corrections to scaling are interpreted in terms of a 
length scale indirectly measuring the knot size [ [Tof . From 
this, it is concluded that knots are only weakly localised, 
in a sense that will be defined below. 

Even if it is computationally quite hard to directly 
quantify the length of a knot in three-dimensional space, 
the same is not true for flat knots which result from the 
projection of a closed curve onto a plane. Such projec- 
tions are well known from knot theory where they are 
at the basis of the determination of topological invari- 
ants jnj. Experimentally flat knots can be realized by 
adsorbing polymers like DNA , or even macroscopic 
chains |pL3|| on a plane. In the former case thermal fluctu- 
ations can allow the system to reach equilibrium in two 
dimensions. 

A flat knot can be drawn as a number of crossings con- 
nected with arcs. At each crossing some monomers have 
to detach from the plane. If the adsorption interaction is 
sufficiently strong, any non-minimal number of crossings 
becomes very unfavourable energetically. For example, 
for the trefoil knot, this minimal number is three and 
there are six arcs. In Ref. jl4| it was shown that for a 
ring of length L with excluded volume and at infinite tem- 
perature, there is typically one of the arcs whose length is 
of order L, whereas the total length of all the other arcs, 
I, is much smaller: / <C L. More precisely, (I), the aver- 
age value of I, was found not to diverge with L. In this 
case one says that the knot is localised. If in contrast, 
(I) ~ L f , one speaks of weak localisation (0 < t < 1) or 
delocalisation [t = 1) of the knot. 

In the present Letter, we investigate the size of flat 
knots when the temperature T is lowered, or equivalently, 
when the quality of the solvent gets worse. Under these 
circumstances the polymer will undergo a collapse tran- 
sition from a coil to a globule shape below a theta point 
temperature Tg Jl5[ . Our main result is that in the col- 
lapsed phase knots are delocalised. At the theta-point, 
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we find them to be weakly localised with t = 3/7. 

A model for flat knots can be defined Jl6| on the square 
lattice whose set of edges is extended with the diagonals 
of the squares. The bonds of the polymer can visit each 
edge and each vertex of this extended lattice at most 
once. A diagonal can only be occupied if at the same 
time the other, perpendicular diagonal within the same 
elementary square is also occupied. Each pair of occu- 
pied diagonals represents a crossing (not to be consid- 
ered as a lattice vertex) in the projection of the knot. 
Thus, one has to further specify which of the two di- 
agonal bonds goes under the other one. The model is 
simulated within the grand canonical ensemble, where a 
fugacity K is assigned to each bond of the ring polymer, 
while the number of crossings is constrained to the min- 
imum consistent with the topology. As usual, to include 
the possibility of theta collapse, we associate an attrac- 
tive energy with each pair of lattice vertices that is visited 
by non-consecutive bonds. Fig. 1 shows a configuration 
with the topology of a trefoil. 




FIG. 1. Typical configuration of a polymer with a trefoil knot 
in the low T collapsed phase. 

In our Monte Carlo approach, a Markov process in 
the configuration space of the polymer is constructed by 
a combination of local and non-local moves. These are 
chosen as in Ref. @] and are such as to ensure invari- 
ance of the polymer topology. Averages at fixed T are 
then calculated using a multiple Markov chain (MMC) 
implementation in the fugacity K jl7]]. This garantees 
an exhaustive sampling of configurations also at rather 
low temperatures First, we obtain precise estimates 
of K C (T), the critical fugacity above which the grand 
canonical average (L) = oo. Since for the case of a flat 
trefoil, there are always only three pairs of diagonals oc- 
cupied, we can expect that Tg is very close to the value 
of an interacting self-avoiding ring model without cross- 
ings. To verify this, we investigated the average squared 
radius of gyration (R^)l f]j| as a function of L for this 
case. Fig. 2 reports our results for different T's. One 
expects that asymptotically 



(-Rg)i, ~ AL 2v (1) 

Our data clearly show the expected three regimes. At 
high T's, we are in a self-avoiding walk regime with v ~ 
3/4. At low T's we determine an exponent v = 0.49 ± 
.02, consistent with the value appropriate for a collapsed 
polymer, v = 1/2. Finally, close to 1/T = 0.67 we find 
a v in agreement with that at the theta-point, i.e. v = 
4/7 HJ. Hence, we estimate 1/Tg = 0.67 ± .02, fully 
consistent with determinations for unknotted rings pi) . 
We also conclude that the exponent v at the theta-point 
and in the collapsed phase is not modified by changing 
the topology from that of a flat unknot to that of a flat 
trefoil. 
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FIG. 2. Log-Log plot of (-Rg)i as a function of L for a ring with 
trefoil knot. From top to bottom : 1/T = 0, 0.50, 0.67(~ 1/Tg), 0.8. 
The dot-dashed, the dashed, and the dotted lines have slopes 3/2, 
8/7, and 1.0, respectively. 

In order to characterise the tightness of the trefoil knot, 
we consider all its six arcs and determine the statistics 
of their lengths lx < I2 < • • ■ < ^6- Clearly the largest 
arc must always have a length proportional to L. In 
Fig. 3, we plot the average length of the second largest 
arc, i.e. (15)l, as a function of L. For high T's, we see 
that (^5)^ saturates at large enough L. Moreover, all the 
other arc lengths remain much smaller, typically only a 
few bonds. We conclude that in the swollen regime the 
knot is localised. The quantity {1)l/L — {{Y^i=\h)) l / L 
approaches zero when L — * 00. 

At Tg, we find instead that grows as T*, with 

t = 0.44 ± 0.02. All other lengths remain again very 
small. Thus, the typical shape of the polymer appears 
to be that of a figure eight, as found also at T — 00 [ pT| . 
We also conclude that the knot is weakly localised. 

Finally, and most interestingly, we find that in the col- 
lapsed phase, and for L-values that are sufficiently large, 
~ L, implying a delocalisation of the knot. In this 
case, there is ample evidence that also the average length 
of smaller arcs, like starts to grow proportional to 
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L for still longer polymer lengths. We suspect that in 
sufficiently long polymers all average arc lengths will be- 
come extensive in L. Thus, a description in terms of a 
figure eight breaks down in the collapsed phase. 
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FIG. 3. Log-Log plot of as a function of L for a trefoil 

knotted polymer. From bottom to top : 1/T = 0,0.50,0.67,0.8. 
The dot-dashed line has a slope 0.44 whereas the dashed line has 
slope 1.0. 

It is possible to associate a polymer network Q with 
N segments 22 1 to each flat knot configuration with the 
same number of non-microscopic arc lengths |H| l . Inter- 
est in polymer networks was revived recently by a suc- 
cesfull development concerning DNA denaturation ]p3[ . 
The same type of approach was subsequently applied in 
fL4| where it was found that flat knots with excluded vol- 
ume at T — oo are always localised. Our data for T > Tg 
support this conclusion, and moreover convincingly show 
that localisation is present throughout the whole high 
temperature phase (Fig. 3). 

On the basis of a network description, this time with 
interacting polymer segments, it is possible to gain fur- 
ther insight into some of our results. Let Q be such a 
network with vertices of degree k connected by N 
arcs of total length L. The partition sum Zg(l\, . . . , In) 
of the network scales as fl^Jzl 



Zg(h, . . . , In) = K c I 
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where Fg is a scaling function and jg = 1 — vdC, + 
Ylk n k (J k- Here d is the dimension of space and C is the 
number of independent loops in the network. The lengths 
of the network segments , corresponding to macroscopic 
knot arcs, are l± < ... < In- Finally, the exponents Cfc are 
anomalous dimensions associated to the fc-leg vertices of 
the field theory describing the polymer in the continuum 
limit p^| . In two dimensions, Coulomb gas methods al- 
low an exact determination of the oVs. In particular, at 
T = oo, where polymers are described by the n —> limit 



of a critical 0(n)-model, one has dfe = (2 — fc)(9fc + 2)/64 
p2|. On the other hand, the theta-point is described 
by the critical low temperature phase of the 0(n = 1)- 
model, for which a k = (2 - k)(2k + l)/42 @. Finally, 
following pa] , we assume that the properties of the low T 
regime can be related to those of dense polymers, which 
are described by the low T phase of the 0(n — 0)-model. 
For this case special care has to be taken, and one ends 
up with a scaling form for Zg which is slightly differ- 
ent from Eq.(H) |£6|. Yet, despite these differences, the 
network picture can still be applied to flat knots with 
cr fe = (4-fc 2 )/32 f§. 

When in a network one or more arcs become very short, 
and hence the crossings on which they are incident ap- 
proach each other very closely, the network should be 
replaced by another, contracted one with fewer segments 
and crossings, and a different jg. In this way, it can be 
understood that even though any projection of a trefoil 
contains six arcs, at a more coarse grained level, the typi- 
cal contribution appearing in a numerical simulation can 
come from a network with fewer crossings. This is what 
happens, for example, at the theta-point, where we found 
that the flat trefoil looks like a figure eight (I4 ~ 0(1)). 
At Tg, where v — 4/7, (0) predicts that the partition 
function of a figure eight network scales as 



Z* =K, 
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with 78 = —12/7. For l/L — » 0, this partition sum should 
in its turn reduce to that of a self-avoiding ring at the 
theta-point, which is known to scale as 2 ~ K~ L L~ vd . 
This simple analysis [g3| then teaches us that F s (x) ~ 
x- c for x -> 0, with c = -(7s - 1 + ud) = 11/7 » 1.57. 
Hence, for l<Lwe predict that 

~ L (L — l)~ vd r c 



k: 



(4) 

In Fig. 4 (left) we present our data forpfis), the probabil- 
ity distribution of Z5, at T = Tg. From (|J) it follows that 
p(h) ~ A fit to this form leads to c = 1.63 ± 0.08, 
consistent with the above prediction. Finally, we get from 
([|) that (Zg) l ~ L 3 / 7 , in good agreement with the numer- 
ical estimate ~ £' 44± 02 . We expect that t — 3/7 is 
an exact result which characterises the weak localisation 
of the knot at the theta-point. 

At T < Tg, an analysis using the results of |2(| can still 
be made for a ring with the shape of a figure eight, and 
leads to the prediction c = 11/8. In Fig. 4 (right), we 
also show our data for p(l 5 ) at 1/T = 0.8 > 1/Tg. There 
is indeed an initial power law decay with an exponent 
1.34 ± .12. But in this case, p(h) flattens for larger Z5- 
values and it is this broadening which eventually leads to 
the delocalisation of the knot. The preasymptotic slope 
indicates that for relatively small I5, the weakly localised 
figure eight network configurations are still dominating 
the partition sum. When (Zg)i ~ L one can get no help 
from network scaling arguments in determining p{l§). In- 
deed, such arguments are only valid for 1$ <Si L which is 
not the relevant range when ~ L. 
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FIG. 4. Log-Log plot of p(l 5 ). On the left (right), 
we show our data at 1/Tg ~ .67 (1/T = .8). The 
dashed lines have slopes of respectively —1.63 and —1.34. 
Different symbols refer to different canonical averages with 
L = 200(A), L = 400(*), L = 600(+), L = 800(D) and 
L = 1400(o). 

We verify that also for other prime knots, such as the 5i 
and 7i |Tl| , the typical network configuration at T = Tg 
is still the number eight and delocalization occurs for 
T < T e . 

In conclusion, we find that a flat prime knot in an ad- 
sorbed polymer remains localised as long as T > Tg, be- 
comes weakly localised at Tg, and eventually delocalises 
at T < T e . While delocalisation is supported by strong 
numerical evidence, other results are fully consistent with 
our predictions based on a polymer network approach. 
One may speculate that the delocalisation found here 
could be a more general phenomenon for the topology 
of interacting polymers. 
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